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We consider the thermodynamics of chiral models in the mean-field approximation and discuss the 
relevance of the (frequently omitted) fermion vacuum loop. Within the chiral quark-meson model 
and its Polyakov loop extended version, we show that the fermion vacuum fluctuations can change 
the order of the phase transition in the chiral limit and strongly influence physical observables. We 
compute the temperature-dependent effective potential and baryon number susceptibilities in these 
models, with and without the vacuum term, and explore the cutoff and the pion mass dependence 
of the susceptibilities. Finally, in the renormalized model the divergent vacuum contribution is 
removed using the dimensional regularization. 
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I. INTRODUCTION 

The critical properties of strongly interacting mat- 
ter near the chiral transition at finite temperature and 
density are studied in first principle calculations using 
lattice simulations of QCD [3|. A complementary ap- 



proach is offered by cfFcctivc chiral models [2|-|llj. like 
the Nambu-Jona-Lasinio (NJL) [3| and the chiral quark- 
meson (QM) [ll[ models as well as the correspond- 
ing extended versions, which include the interaction of 
quarks with a uniform temporal gluon field. The lat- 
ter, the Polyakov loop extended Nambu-Jona-Lasinio 
(PNJL) d g and quark-meson (PQM) 0, Q models 
reproduce prominent features of QCD thermodynamics, 
obtained in lattice simulations. In particular, both mod- 
els exhibit a chiral as well as a dcconfinement transition. 
In this paper, the consistency of mean-field calcula- 
tions of many-body systems within chiral effective models 
like the NJL/PNJL and QM/PQM models is addressed. 
In particular, wc discuss the role of the lowest order 
fermion vacuum contribution. In the NJL/PNJL model, 
this term is responsible for the dynamical breaking of the 
chiral symmetry in vacuum and is therefore necessarily 
taken into account. On the other hand, in the QM/PQM 
model, the spontaneous breaking of the chiral symmetry 
is due to the meson potential. Hence, for qualitative con- 
siderations, the fermion vacuum loop is frequently omit- 
ted p, [11413 • 

We show that a mean-field approximation, where the 
fermion vacuum term is neglected, commonly referred to 
as the no-sea approximation, leads to a distortion of the 
critical behavior at the chiral transition. In particular, 
in this approximation the order of the transition in the 
chiral limit is always first order. In contrast, when this 
term is included, the transition is of first or second order, 
depending on the choice of coupling constants and on the 
baryon density. 

Differences in the thermodynamics between the NJL 



and QM models found in the literature can be attributed 
to the use of the no-sea approximation in the QM model. 
For instance, in this approximation, the adiabatic tra- 
jectories obtained in the QM model exhibit a kink at 
the chiral crossover transitiori, while in the NJL model 
they are smooth everywhere [1^. This effect is a trace 
of the underlying first-order phase transition of the QM 
model in the chiral limit, when fermionic vacuum fluc- 
tuations are neglected ^M,- Moreover, the omission of 
fermion vacuum fluctuations can result in a violation of 
the semi-positivity of the spectral function as well as in 
a violation of detailed balance [13 • The vacuum fermion 
fluctuations were also shown in the PQM model to con- 
siderably affect the structure of the phase diagram in the 
presence of an external magnetic fleld [201 ■ 

In this paper we demonstrate the importance of 
fermion vacuum fluctuations for a consistent formula- 
tion of the thermodynamics of the QM and PQM models 
in the mean-fleld approximation. We show that, in or- 
der to reproduce the second-order chiral phase transition 
expected for two flavor QCD within these models j2]| . 
it is necessary to take the fermion vacuum term in the 
thermodynamic potential into account and to perform 
a proper regularization of the divergence. Furthermore, 
we also show that when the chiral symmetry is explicitly 
broken, the thermodynamics is strongly affected by the 
no-sea approximation. Finally, we remove the inflnities 
of the vacuum term using the dimensional regularization. 



II. THERMODYNAMICS OF THE 
QUARK-MESON MODEL 

The chiral QM model is used as an effective realiza- 
tion of the low-energy sector of QCD. The replacement 
of the SU{Nc) local gauge invariance of QCD by a global 
symmetry in the QM model eliminates the possibility to 
address conflnement within this model. However, as re- 



ccntly shown, the confining properties of QCD can be ef- 
fectively accounted for by introducing the Polyakov loop 
in the chiral quark-meson model [3] , in complete analogy 
to the PNJL model p. Consequently, the PQM model is 
an effective model of both the chiral and confining prop- 
erties of QCD. 

The Lagrangian of the PQM model reads 



'U {a, tt)-U {£,£*) 



(1) 



where U{£,£*) is an effective potential of the gluon field 
expressed in terms of the thermal expectation values of 
the color trace of the Polyakov loop and its conjugate. 



with 



e = ±{TrMS)), r = ^(TreLt(f)), (2) 



L{x) — V cxp i / dTA4{x,T) 



(3) 



Here V stands for the path ordering, (3 — 1/T and A4 ~ 
iAq. We use the short-hand notation A^, = g\°-A''^/2. 

The 0(4) representation of the meson fields is (/> = 
(cr, tt) and the corresponding SU{2)l (^ SU{2)ji chiral 
representation is given by cr -f i^^f ■ if. This implies that 



J 



there are iV? = 4 mesonic degrees of freedom coupled to 
Nf = 2 flavors of constituent quarks. The coupling be- 
tween the effective gluon field and quarks is implemented 
through the covariant derivative 



Df, = df,, - iAf,, 



(4) 



with Afj_ = SfioAo. The purely mesonic potential of the 
model reads 



A 
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U{a,Tf)='-{a^ + TT^^v^y -ha, (5) 



while the Polyakov loop potential is parameterized by a 
Z{3) invariant form: 



UjT-JJ*) _ b2{T) 
T4 2 

The parameters. 



ti- 



-^{f'+t^) + ^{ttf . (6) 
6 4 



b,{T) = a, + a,[^\+aJ^\\aJ^\\ (7) 



where a^ — 6.75, ai = —1.95, 02 ~ 2.625, 03 = —TAA, 
63 — 0.75 and 64 — 7.5 are chosen to reproduce the equa- 
tion of state obtained on the lattice for pure SU(3) gauge 
theory [6|, which yields a first-order phase transition at 
the temperature To = 270 MeV. 



The thermodynamic potential of the PQM model is in the mean-field approximation (7[ 



nMFiT,fi;{a),£,n^^iT;e,n + U{{<j),{7T)=0) + ngg{T,fi;{a),£,n 
The quark contribution, with dynamical mass rriq = g{(T 



f],,(T,/.; {a),£,£*) = -2Nf / —^ {N^E.eiA^ - f) 



(8) 
and energy Eg — \ p"^ -t- m^, is given by 

+ rin.gW(T,/i;(a>,£,r) + rin.g(-)(r,M;(a),£,r)}, (9) 



where 



g^+\T,^i-{a),£,t) = l + Uer^^^' 
g^~\T,^l;{a),£,£*) = g^+\T,-fi; 



a), £*,£). 
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(10) 
(11) 



The first term in Eq. © is the fermion vacuum contri- 
bution, regularized by the ultraviolet cutoff A. 

The equations of motion for the mean fields (a) , £ and 
£* follow from the stationarity conditions 

oUmF diiMF c/sjmf 



= 0, 



(12) 



da d£ d£* 

which yield the temperature and chemical potential de- 
pendence of the chiral (cr) (T, /x) as well as of the Polyakov 
loops £{T, ^) and £*{T, ^) order parameters. 



The parameters h, g, A and v are specified by requiring 
that the following vacuum properties be reproduced: the 
chiral condensate is fixed to (c)(0, 0) ~ /^r, the pion mass 
to mj = h/ f-ji, the constituent mass of quarks to rrig = 
gf^: and the sigma mass to 



2 a2i7MF(o,o;/,,^,r) 



9fj2 



(13) 



at the potential minimum. In the following we use the 



700 McV, m° 



335 MeV 



values mjr = 138 McV, rua 
and /^ = 93 MeV. 

The results in the chiral limit arc obtained with the 
same parameters except for the symmetry breaking term 
h, which is set to zero. Thus, we ignore any influence of 
the pion mass on the other parameters, like e.g. f-^ [22 |. 
For the problem at hand, such shifts are unimportant. 

The thermodynamic potential in the QM model is ob- 
tained from Eq. ^ by setting £ — £* — ! and ignoring 
the constant contribution from the Polyakov loop poten- 
tial. In this case the functions g^^'' reduce to 



W(T,/x; (a),£ = l,r = 1) = [1 - nFiT,±ii; (a))] 



9 



where 



nF{T,±n;{a)) = 



exp I ^^^ I + 1 



(14) 



is the Fermi-Dirac distribution function. In the following, 
we do not consider fluctuations of the a field. Hence, we 
simplify the notation by dropping the bra-ket notation. 
Thus, from this point on, a = (a) denotes a uniform 
classical field. 



III. LANDAU EFFECTIVE THEORY 



where 'Je is Euler's constant. Note that there is no cu- 
bic quark mass term in this expansion in contrast to the 
bosonic case. 

The effective potential of the QM model is the sum of 
the quark and meson contributions 



n{a) = n 



th 



\V^ 



A 



(16) 



We note that for A > 0, one would expect the phase 
transition to be second order. 

Because the dynamical quark mass is proportional to 
the order parameter a, the logarithmic term in Eq. ()15|) 
competes with the quartic term in the meson potential 
P^ . For sufficiently small values of the order parame- 
ter, the logarithmic term is arbitrarily large and negative. 
This leads to a change in sign of the effective quartic cou- 
pling in Eq. p6p . which implies that the phase transi- 
tion is first order for any finite value of A. This effect re- 
sembles the Coleman- Weinberg fluctuation-induced first- 
order phase transition [2g at T = 0. In order to make the 
mechanism more transparent, we introduce the Landau 
effective potential, corresponding to Eq. (|16p 
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n{a) - 17bg = -a{T)a' + -ba^[l + cln(a/ao)], (17) 



The mean-field thermodynamic potential ([8]) can be 
used to explore the infiuence of fermion vacuum fiuctua- 
tions on the critical and thermal properties of a system of 
quarks interacting with the meson and gluon mean fields. 
In this section we explore the role of the fermion vacuum 
term in an effective Landau theory for the QM model 
with and without the fermion vacuum term. Thus, for 
the moment we neglect the Polyakov loop potential and 
its coupling to the quarks. The conclusions of this sec- 
tion are, however, directly applicable also to the Polyakov 
loop extension of the QM model. 



Effective potential without fermion vacuum 
fluctuations 



The contribution of the fermion vacuum term to the 
effective potential ^ can be eliminated by setting the 
cutoff parameter A to zero. In the chiral limit, i.e., for 
/i — >■ 0, and near the chiral phase transition the dynamical 
quark mass, which is proportional to the order parame- 
ter, is very small. Consequently, the thermodynamics 
of the QM model can be formulated as an effective the- 
ory for the order parameter. Due to the smallness of 
the quark mass the high temperature expansion can be 
applied to derive the effective quark potential 1231 - 1251 



n 



th 



NcNfT 



.n(5)... 




where, to linear order in T — Tc, a{T) = A- {T — Tc) and 
A> 0. The quartic coupling b as well as the coefficient c 
are assumed to be positive constants. In general, both b 
and c depend on temperature. However, this dependence 
is irrelevant for our discussion, as long as b{T(^ and c{T(^ 
are non-zero. 

The parameters of the Landau potential (flT)) can be re- 
lated to that of the QM model. By inspection of Eq. (fT5|) 
and (1161) one finds 



6 



(T' " T^) 



A-{T~Tc) + {T -T, 



where 



T: 



6Aw2 



g^NcNf ' 
A = g^NcNfTc/3. 



(18) 

(19) 
(20) 



Similarly, one can identify the quartic coupling b and the 
coefficient c 



6 = A 



g^NJVf_ 
27r2 



lE-^ 



2TT^b ' 



(21) 
(22) 



as well as the parameter ctq i which for T ^ Tc is given by 



c^o 



9 



a 



(23) 
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FIG. 1: The effective potential as a function of tfie order parameter in the vicinity of the pliase transition for vanishing (left 
panel) and finite (right panel) cutoffs. The curves correspond to /i = and a few temperature close to TpT- The potential is 
normalized to zero at the origin. 



The background contribution llbg in Eq. \YJ\) is indepen- 
dent of the order parameter a^ and given by the T^ term 
in Eq. (dSl). 

For small positive values of a, i.e. for T > T^-, the 
potential (|17[) has two minima located at 



CTl/cTo = 0, 

o-2/fTo = K 



2bcKal ' 



(24) 
(25) 



where k = exp(— 1/4 — 1/c). At these points, the corre- 
sponding values of the potential are 






(26) 
bc{aoK)'^ + a{aoKf. (27) 



Consequently, for a < jbc{<jQK)'^ the second solution, 
tT2, is energetically favored. For a = jbc{aoK)^ the two 
minima at cti = and CT2 = f'^o'o are degenerate, i.e. 
fl{ai) ~ ^(cr-i) ~ f^bg- Finally, for a > jbc{<7oK)^ the 
first solution ai is energetically favored. This analy- 
sis indicates that a first-order phase transition occurs at 
TpT = Tc + &c(k(To)^/(4A), i.e. at a temperature above 
the nominal critical temperature Tc. 

The results obtained above are illustrated in Fig. [!]■ left, 
where we show the effective potential of the QM model 
for /i = 0, neglecting the fermion vacuum term. The de- 
pendence of the potential on the sigma field, for selected 
values of the reduced temperature t = {T — Tpt)/Tpt, 
shows explicitly that this model exhibits a first-order 
phase transition in the chiral limit |33| . Thus, without 
the fermion vacuum term, this model does not reproduce 
the second-order transition of QCD with two massless 
fiavors, expected from the universality argument. 



B. Effective potential with the fermion vacuum 
contribution 



We now include the contribution of the fermion vac- 
uum term of the QM model and analyze the resulting 
thermodynamic potential. For sufficiently small values 
of the dynamical quark mass niq, the vacuum term in 
Eq. ^ can be expanded in a Puiscux series in mq/A. 
The leading terms of this expansion reads 



99 



Eq9{A^-p^) 




(28) 



4 In 



2A/ 
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The logarithmic dependence on niq cancels between the 
thermal and vacuum contributions to the thermodynamic 
potential, ([T51) and ([25)1 . Consequently, when the fermion 
vacuum term is included, the effective potential has the 
typical structure expected for a theory with the second- 
order phase transition 



Q(a) - Ob 



1 



a{T)a^ 



1 



ba\ 



(29) 



The mass coefficient is formally given by the Eq. (|18p . 
however with the following redefinition of the critical 
temperature 



(30) 



^ 7r2 


A 2 6Aw2 

g^NcNf 


the quart ic coupling reads 


--"tf 


. , /^2A\1 

7_E - 1 + In -— 
K'^TcJ 



(31) 



In the above expressions the parameters A and v'^ de- 
pend on the cutoff A, because we require that the physi- 
cal vacuum parameters discussed above should be repro- 
duced independently of the cutoff. Taking this depen- 
dence into account, the leading A^ term in Eq. ((30)) is 
canceled by a corresponding term in Xv^. Consequently, 
for large A, the critical temperature is cutoff indepen- 
dent: 



r2(A > Tc 






3{r 



27r2 



+ 0{A-^). (32) 



The quartic coupling is not positive definite for all 
values of A. For A > A*, where A* is defined by 
b{T -^ Tc,A*) = 0, the thermodynamic potential ((29)) 
has only one minimum at any temperature, at cti = or 
at (72 = a/— a/6. Thus, a second-order phase transition 
takes place at TpT — Tc- This behavior is illustrated 
in Fig. [l]-right, where we show the thermodynamic po- 
tential obtained in the QM model, including the fermion 
vacuum term, for A = 0.6 GeV, which implies A = 3.06 
and v^ = —(578)2 MeV^. For this choice of parameters, 
the position of the minimum of the thermodynamic po- 
tential varies continuously with temperature, indicating 
a second-order phase transition. 

However, for A < A*, the quartic coupling is negative, 
while, as shown in Eqs. ((15)) and ((28)) . the coefficient of 
the a^ term is always positive, as required for stability. 
The negative value of b results in the appearance of the 
second minimum in the potential, leading to the first- 
order phase transition. This is similar to the situation 
described in the previous section, where the fermion vac- 
uum term was neglected. The appearance of a first-order 
phase transition for small A is illustrated in Fig. [21 where 
the temperature dependence of the order parameter is 
shown for several values of the cutoff A. For A < 0.3 
GeV, the order parameter exhibits hysteresis, character- 
istic of a first-order phase transition. For the parameter 
set employed here, we find A* ~ 244 MeV. This value 
is comparable to the physical scales that characterizes 
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FIG. 2: The order parameter as a function of temperature 
for different values of the cutoff parameters. The uppermost 
line corresponds to the renormalized theory. 



the vacuum. In general, the cutoff should be larger than 
any physical scale of a system. In this case we find a 
second-order chiral phase transition in the QM model. 

The fermion vacuum fluctuations modify the critical 
properties of the QM model considerably. The modifi- 
cation depends on the value of the ultraviolet cutoff A, 
which is introduced to regularize the divergent vacuum 
term. For large cutoffs (A > A*), the QM model ex- 
hibits a second-order phase transition, whereas in the 
opposite limit the phase transition stays first order, as in 
the approximation, where the vacuum term is completely 
neglected. 

Clearly, a finite cutoff modifies not only the critical 
structure of the medium but also the properties of physi- 
cal observables. Hence, it is desirable to remove the cut- 
off dependence by means of a suitable renormalization 
procedure. This problem will be addressed in the next 
section. 



RENORMALIZATION OF THE FERMION 
VACUUM TERM 



IV. 



The expansion ((28)) indicates that, in addition to the 
fermion vacuum term, also the mass parameter a and 
the quartic coupling h diverge in the limit of A — >■ oo. 
These divergencies can be regularized by introducing a 
cutoff, as done above. In this section we choose a more 
convenient route and use the dimensional regularization 
scheme. 

The vacuum term is, to lowest order, just the one-loop 
effective potential at zero temperature 24 1 



Qvac 
"99 



'2NfN, 



d^p 



3^9 



= -2NfN,j ^^HpI + ED 



+ C, (33) 



where the infinite constant C is dropped, since it is inde- 
pendent of the fermion mass. 

To regularize Eq. p3p we perform the dimensional reg- 
ularization near three dimensions, d — i ~ 2e. The re- 
sulting potential up to zeroth order in e reads 






1 



-3 + 27£ 



4 In 






-^ 



\2^M J 

(34) 

where M is an arbitrary renormalization scale parameter. 

The thermodynamic potential is then renormalized by 

adding a counter term to the Lagrangian of the QM or 

PQM model. With the choice [sj] 

^^ = ^.9^-'{7-^h3 + 2,,-41n(2V¥)]}, 

(35) 
the renormalized contribution of the fermion vacuum 
loop reads 



"«« - 87r2 ™«^''U/ 



(36) 



The quark contribution to the thermodynamic potential 
is then given by Eq. Q, with the vacuum term replaced 
by Eq. Ipj). 

For large values of the order parameter, nig = ga > M , 
the renormalized vacuum contribution in Eq. p6p is un- 
bounded from below. At T = 0, the potential ([5|), includ- 
ing the renormalized quark vacuum term, has a minimum 
at the physical point cr = /tt, a maximum at a ~ 212 
MeV and then decreases monotonously to — oo for cr — > oo 
[SSt - The instability for large a is symptomatic of the 
renormalized one-loop approximation j2a|.The inclusion 
of higher order loop contributions is known to cure this 
problem. 

Dimensional regularization introduces an arbitrary 
renormalization scale parameter M. Consequently, the 
parameters v and A are M dependent through the condi- 
tions (|13p and p2)) . Nevertheless, in the one- loop approx- 
imation, physical quantities are independent of M [26j . 
Indeed, it is rather straightforward to show that 



dnMF{^[M],v^[M],M) 
dM 



(37) 



SufBciently close to the phase transition, the Landau 
potential corresponding to the renormalized QM model 
has the form (|29p . with the mass parameter given by 
Eqs. (|18ll20p . The critical temperature T^ is independent 
of the renormalization scale, since d{Xv'^)/dM = 0. Con- 
sequently, it can be expressed solely in terms of physical 
quantities. 



rri2 



3m2/2 



NcNfimOr 



3(m^ 
27r2 



(38) 



This expression for Tc agrees with the one obtained in 
the cutoff regularization scheme, Eq. (|32p . for A — > oo. 
Also the quartic coupling b can be expressed in terms of 
physical quantities 



6 = 






27r2 



7-E 



In 



■kTc 



(39) 



Thus, we conclude that in the mean-field approxima- 
tion, the thermodynamic potential of the quark-meson 
model, including the renormalized fermionic vacuum con- 
tribution, is uniquely defined. In Fig. [5] we show the 
temperature dependence of the order parameter, as ob- 
tained within this model. This and other physical quanti- 
ties are, as indicated above, free from ambiguities related 
to the choice of the renormalization scale. In the chiral 
limit, this model exhibits a second-order phase transi- 
tion at /i = 0, as expected for two-flavor QCD. Thus, 
the quark-meson model, with appropriately renormalized 
fermionic vacuum fluctuations, is an effective QCD-like 
model, suitable for describing the physics at the chiral 
phase transition. 
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FIG. 3: The kurtosis in the PQM model as a function of 
temperature for several values of cutoff parameters A. 



V. QUARK NUMBER DENSITY 
FLUCTUATIONS 

As discussed above, the fermion vacuum term is re- 
quired to obtain the correct critical behavior in the chiral 
limit. In this section, we show that the vacuum term is 
important also for the thermodynamics for non-zero pion 
masses. To illustrate this, we consider fluctuations of the 
net quark number density within the PQM model. 

The fluctuations of the quark number density are char- 
acterized by the generalized susceptibilities. 



.m- 



a"b(T,A^)/r4 



difi/TY 



(40) 



M=0 



The first two coeffients, C2 and C4, are the second and the 
fourth cumulants of the net quark number Nq, 



C2 



C4 



y2 



{m,r 



{{6Ngr)~3{{SNg)Y. 



(41) 
(42) 



Here Xg is the quark number susceptibility and SNq = 

Ng^iNq). 

The coefficients c„ (T) are sensitive probes of the chiral 
phase transition, see e.g. [29|. In particular, C2 diverges 
at the critical endpoint of the QCD phase diagram [2[. 
In the chiral limit and at non-zero chemical potential, all 
generalized susceptibilities Cn(T) with n > 2 diverge at 
the 0(4) chiral critical line [30[. Moreover, if the spinodal 
instability is reached, fluctuations of the net quark num- 
ber also diverge at a first-order chiral phase transition 

i. 

A very particular role is attributed to the so called 
kurtosis of the net quark number fluctuations [3, [2^, |30, 
l3ll | , the ratio of C4 and C2 



^4.2 — — ■ 



14,2 



C2 



(43) 



This observable is not only sensitive to the chi- 
ral phase transition but also probes the confinement- 
deconfinement transition. At very low and very high 
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FIG. 4: The second C2 and the fourth a order cumulants, computed in the PQM model including the renormalized fermion 
vacuum potential, as functions of temperature for several values of the pion mass. 



temperatures, the kurtosis reflects the quark content of 
the baryon-number carrying effective degrees of freedom 
[30l |3l|. Thus, at low temperatures, in the confined 
phase, i?4.2 — N'^ = 9 while in an ideal gas of quarks [3g 
Ri,2 '^ 1- The influence of the chiral phase transition on 
the kurtosis has been studied within lattice gauge the- 
ory [32]. Due to the "statistical confinement" realized in 
the PQM and PNJL models, these models can also be 
used to explore the interplay between confinement and 
chiral restoration e.g. in the kurtosis |8l.[3]|. 

We first consider the dependence of the ratio i?4_2 on 
the cutoff A in the PQM model, including the vacuum 
contribution, Eq. (P5|) . In Fig. [3] the kurtosis is shown 
at the physical pion mass, as a function of the temper- 
ature for several values of A. The kurtosis converges to 
Ri.2 = 9 at low temperatures and to i?4^2 '^ 1 in the 
high temperature limit for all values of cutoff parame- 
ter. However, near Tc the kurtosis is very sensitive to 
A; for small cutoffs it develops a peak, which is maximal 
for A — > 0. The sharp peak, obtained for vanishing cut- 
off, is a consequence of the underlying first-order phase 
transition in the chiral limit. Clearly, such a strong de- 
pendence on the ultraviolet cutoff is unphysical and must 
be removed by a suitable renormalization of the vacuum 
term. 

Having established the relevance of the fermion vac- 
uum term for the kurtosis, we now study fluctuations of 
the net quark density in the PQM model, including the 
renormalized vacuum contribution. In Fig. 2] we show 
the second and fourth cumulants, C2 and C4, as functions 
of temperature for several values of the pion mass. For 
small 7n^, C2 exhibits a kink at the chiral crossover tran- 
sition, while C4 shows an rapid drop, which in the chiral 
limit evolves into a discontinuity. 

These results can be understood in Landau theory. For 
small a^ the thermodynamic potential reads 



1 



1. 



where flbg{T, fi) is the background (regular) contribution, 
which is independent of a. The external field h breaks 
the symmetry and generates a finite pion mass. For a 
small /i and T near the chiral transition temperature, 
the coefficient a(r, ^i) can be parameterized by [3^ 



aiT,fi)=A 



T-Tc 



-'1 



(45) 



where both, A and B are positive constants. The coeffi- 
cients of the effective potential can be expressed in terms 
of the parameters of the (P)QM model by considering the 
expansion of the thermodynamic potential from Eq. ([8]) 
for small m^ and a. 

In the chiral limit and at yu = the singular contribu- 
tion to the second and fourth cumulants is 



cf^ - ^^^e(r.-T) 



bT^ Tr 



6^2 



0{Te, - T) 



(46) 
(47) 



n{T, /i) - a,g{T, /i) = -a{T, fi)a' + -ba^ - ha, (44) 



Hence, C2 has a kink at the critical point, while C4 is 
discontinuous, in agreement with the numerical results 
for the PQM model presented above. 

The peak structure in C4 near the chiral transition, 
shown in Fig. 21 results from an interplay of the regular 
part, which increases as the critical point is approached, 
and the discontinuity of the singular part. 

The dependence of the fourth cumulant C4 on the pion 
mass is also reflected in the kurtosis of the quark number 
fluctuations. In Fig. [5] we show i?4.2 in the vicinity of 
the chiral phase transition, for several values of the pion 
mass. The transition between the low and high temper- 
ature regimes is sharpened as the pion mass is lowered, 
indicating that this observable is sensitive not only to 
deconflnement but also to chiral dynamics. In the mean- 
field approximation, a clear enhancement beyond the low 
temperature value, i?4.2 = 9, is found just below the 
pseudo-critical temperature, at the physical value of the 



pion mass. There is a further increase of the enhancement 
when approaching the chiral hmit, where i?4_2 exhibits 
a discontinuity at the chiral phase transition, in agree- 
ment with Eq. (j47p . However, when the renormahzed 
vacuum term is included, this dependence is weaker than 
previously obtained within the PQM model, neglecting 
vacuum fluctuations |3ll |. This shows, that an appropri- 
ate rcnormalization of the fermionic vacuum potential is 
essential for a consistent formulation of the thermody- 
namics of the PQM model. 



VI. SUMMARY AND CONCLUSIONS 

We have discussed the importance of the divergent 
fermion vacuum potential in a consistent formulation of 
the thermodynamics of the quark-meson (QM) and the 
Polyakov loop extended quark-meson (PQM) models. 

We have shown that the vacuum term influences the 
order of the chiral phase transition in the chiral limit. 
The QM and PQM models exhibit a first-order phase 
transition at /i = 0, when the fermion vacuum term is 
neglected. On the other hand, when fermion vacuum 
fluctuations are included, a second-order phase transition 
of the 0(4) universality class is possible, as expected in 
models that exhibit the SU{2) x SU{2) chiral symmetry. 

Also for the physical value of the pion mass, the vac- 
uum term influences thermodynamic observables near 
the chiral crossover transition. We have illustrated this 
explicitly, by computing the second and the fourth cu- 
mulants as well as the kurtosis of the net quark num- 
ber fluctuations for several values of the pion mass using 
different regularization schemes for the fermion vacuum 
potential. 



Finally, we have renormahzed the fermion vacuum 
term to obtain the thermodynamic potentials of the QM 
and PQM models, free of any dependence on the unphys- 
ical ultraviolet cutoff. 
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FIG. 5: The kurtosis _R4,2 as a function of temperature for 
several values of the pion mass in the PQM model, including 
the renormahzed fermion vacuum term. 
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